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Abstract 

^ We discuss the physics of vortices in the circumsteUar disks associated with young 

steUar objects. We elucidate the basic physical properties of these localized storm sys- 
tems. In particular, we consider point vortices, linear vortices, the effects of self-gravity, 
, ^ magnetic fields, and nonlinear aspects of the problem. We find that these vortices can 

^ exist in many different forms in the disks of young stellar objects and may play a role in 

y—i the formation of binary companions and/or giant planets. Vortices may enhance giant 

planet formation via gravitational instability by allowing dust grains (heavy elements) 
to settle to the center on a short timescale; the gravitational instability itself is also 
enhanced because the vortices also create a larger local surface density in the disk. In ad- 
dition, vortices can enhance energy dissipation in disks and thereby affect disk accretion. 
Finally, we consider the possibility that vortices of this type exist in molecular clouds and 
On in the disk of the galaxy itself. On all of these size scales, vortices can produce long-lived 

structures which may correspond to observed structures in these systems. 
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|-j Subject headings: stars: formation - hydrodynamics - planets: formation 
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1. INTRODUCTION 

CircumsteUar disks play an important role in the star formation process (see, e.g., 
the reviews of Shu, Adams, & Lizano 1987; Bertout 1989; Beckwith & Sargent 1993). 
Through many recent studies, we now have a reasonably good understanding of the basic 
physical properties of these disks, such as total disk mass, temperature distributions, and 
radial size (see, e.g., Rydgren & Zak 1985; Rucinski 1987; Kenyon & Hartmann 1987; 
Adams, Lada, & Shu 1987, 1988; Beckwith et al. 1990; Adams, Emerson, & Fuller 1990). 
Unfortunately, however, many different physical processes occur in these systems and it 
is not yet known which processes are most important for disk evolution (see, e.g., the 
review of Adams & Lin 1993). In this paper, we explore a new type of physical process 
in circumsteUar disks. We study the basic physics of vortices, which are essentially long- 
lived storm systems. We explore possible effects that vortices can have on disk evolution 
and dynamics. In particular, we show that these vortices can enhance the formation of 
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giant planets within the mechanism of gravitational instability. In addition, vortices can 
enhance energy dissipation and hence affect disk accretion. 

The most well known and well studied example of a fluid vortex in an astrophysical 
setting is the Great Red Spot of Jupiter (e.g., IngersoU 1990; see also the special issue of 
Chaos 1994). This stable swirling storm system has lived for many dynamical timescales. 
Although the exact nature of the Great Red Spot remains somewhat controversial (see, 
e.g., Marcus 1993 and Petviashvili & Pokhotelov 1992 for different perspectives), the 
basic physics is known and provides a starting point for this present work. In particular, 
in this paper we shall focus on quasi-geostrophic vortices, leaving more complex vortex 
models for future work. Similar vortices arise in most rotating fluid systems, such as 
the atmosphere and oceans of the Earth (e.g., Hess 1959; Ghil & Childress 1987). We 
note that vortices have also been studied in the context of galactic disks (Korchagin 
& Petviashvili 1985, 1991; Korchagin, Petviashvili, & Ryabtsev 1988). These previous 
studies have not included the self-gravity of the system and have not considered nearly 
Keplerian rotation curves such as those found in circumstellar disks. We note, however, 
that vortices in circumstellar disks have been discussed previously in a less modern context 
by Von Weizsacker (1944). In addition, two recent papers (Tanga et al. 1995; Barge & 
Sommeria 1995) also discuss planet formation and dust settling in vortices. 

This paper is organized as follows. We flrst formulate the general problem of a 
vortex in a self-gravitating circumstellar disk (§2). Our goal is to flnd analytic (and semi- 
analytic) solutions for the vortices. We consider several different limiting cases in order 
to elucidate the basic underlying physics of these storm systems. We consider the case of 
point vortices in §3; we also present results for linear vortices without self-gravity. Next, 
we consider higher order effects including self-gravity, nonlinear effects, and magnetic 
flelds (in §4). We discuss applications of vortices to astrophysical problems in §5. We 
conclude, in §6, with a summary and discussion of our results. 



2. GENERAL FORMULATION 

In this section, we derive the basic equations which describe the physics of self- 
gravitating vortices in circumstellar disks. In order to obtain analytic results, we make a 
number of approximations along the way. We consider the vortices to be small compared 
to the radial position r in the disk; we thus invoke a local approximation. We can 
then separate the vortical motion from the mean flow of the disk. We introduce both 
the vorticity and the stream function and obtain equations in terms of these physical 
quantities. 

We begin by writing the fluid equations in a form appropriate for a circumstellar 
disk. The radial force equation becomes 

du du V du dh dtb 

^+^^^ + -^-- + ^ + ^ = 0, 2.1 

Ot or r ocp r or Or 

where u {v) is the radial (azimuthal) velocity, ij) is the gravitational potential, and h is 
the enthalpy. We have thus assumed a barotropic equation of state for the fluid (see. 
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however, §4.5). The azimuthal force equation is 



c?u ^ c?u ^ u c?u ^ ^ 1 5 ^ ^\ ^ (2 2) 

dt dr r d(f) r r d(f>\ J 

The continuity equation takes the form 

da 1 d , , 1 d , , , , 

— + -^{rau) + -^{(7v) = 0. 2.3 
Ot r or r ocp 

The gravitational potential ij) obeys the usual Poisson equation, i.e., 

VV = 47rGp. (2.4) 

2.1 Local Approximation 

We want to consider disturbances which are small compared to the radial extent of 
the disk. As a general trend, the disturbance size will be roughly comparable to the scale 
height H of the disk and for thin disks if/r ^ 1. We thus work in terms of the variables 
X and y centered on a point (ro, ^o) rotating with the disk: 

r = ro + X , (2-5) 

y = r((^-(^o) = r((^-Oot), (2.6) 
where Oq is the rotation rate at the center of the vortex (see Figure 1). 
In terms of these new variables, the equations of motion become 

ut - rfloUy + uuj, + vuy - — + + /ij. = , (2.7) 

uv 

vt - rQoVy + uvx + vvy -\- — -\- tpy -\- hy = , (2.8) 

at - rfloCTy + dj,{ua) + dy{va) = . (2.9) 

2.2 Mean Flow Approximation 

As the next simplification, we will separate out the mean flow from the disturbance. 
In equilibrium, the disk is azimuthally symmetric and has no radial flow, i.e., u = and 
V = vq = rO(r). Here, we want to separate the perturbed flow of the vortex from the 
unperturbed motion; we therefore write the azimuthal velocity v as the sum of parts: 

V = vo + vi = rQ{r) + vi , (2.10) 

where 0(r) is the unperturbed rotational speed as a function of radius. This separation 
of the mean flow can be done in a variety of ways and can lead to somewhat different- 
looking equations of motion (see, e.g., Petviashvili & Pokhotelov 1992; Ghil & Childress 
1987). 
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Next we define an operator I), 

V = dt + [v^- rOo]a, = dt + [rO(r) - rOo]a, . (2.11) 
The equations of motion now become 



T>u + uux + viUy — 2Vtvi -\- hij. -\- tpix = , (2.12) 

r 

Vvi + uvix + viviy + uO + — !- + uvox + hiy + tpiy = , (2.13) 

Va + d^{ua) + dy{via) = . (2.14) 

In obtaining this form for the force equations, we have used the fact that the mean 
(unperturbed) flow is determined by the equation 

^ = rQ^ = d^ibo + ho) . (2.15) 
r 

We have thus separated out the unperturbed part ho and the remaining part hi of 
the enthalpy function. We have also invoked this same separation for the gravitational 
potential, i.e., ?/> = ^/jq + V'l- Although we have separated out the perturbed quantities 
from the background, notice that we have not assumed that the perturbations are small. 
Thus far, we have retained the full nonlinearity of the equations of motion. 

Next, we can combine the two components of the force equation into a single vector 
equation. Keeping only the leading order terms in x/r, we thus obtain 

Dv + (v • V)v + 20(5 X v) + V/ii + VV'i + ur(dM)y = , (2.16) 

where v = (u, ui) is the perturbed velocity fleld. Notice the presence of the Coriolis term 
20(5 X v), which arises because we are working in a rotating frame of reference. 

We also want to write the continuity equation in terms of enthalpy instead of in 
terms of the surface density a. The enthalpy is deflned by 

V/i=— V(7, (2.17) 

(7 

where a is the effective sound speed of the fluid. We can thus write the continuity equation 
in the form 

D/ii +v- V/i + a^V- V = 0, (2.18) 

where we have used the fact that the unperturbed state is independent of both time t 
and the variable y so that T>h = T>hi . 

2.3 Vorticity 

Next, we want to write the equations of motion in terms of vorticity, instead of in 
terms of velocity. This approach makes sense because we are interested in vortex solutions 
and vorticity is the relevant physical quantity for vortices. The vorticity io is deflned by 

cu = 5-(Vxv), (2.19) 



4 



i.e., we consider only the z component of the vorticity. 

We now take the curl of the force equation (in the form of equation [2.16]). Next, we 
combine the continuity equation [2.18] with this newly derived vorticity equation. After 
some rearrangement, this equation of motion becomes 



D 



+ 20 + rd^9. 



a 



(2.20) 



where we have defined the operator D to be the total advective derivative operator: 



D = D + V • V = a* + r(0 - no)dy + v • V 



(2.21) 



2.4 Geostrophic Balance 

We next use the fact that the fluid can be in some kind of generalized geostrophic 
balance. Geostrophic balance occurs when pressure forces are balanced by Coriolis forces 
(so that the inertial forces are negligible). In this work, we generalize this concept to 
include gravity, so that both pressure and gravitational forces are balanced by the Coriolis 
force (see Figure 2). In particular, we write the vector velocity v in the form 

^ [5 X (V/ii + VV-i)] , (2.22) 



20„ 

where Oq is the value of the rotation rate in the disk at the location of the origin of the 
local (x,y) coordinate system. In the geostrophic approximation, the vorticity can be 
written 

io = -^{V^h+V\h)- (2.23) 
Thus, our equation of motion [2.20] now becomes 



D 



^ .{V'h,+A7rGp^) + ^^^±^]=0, (2.24) 



2O0C c 



where we have used the flrst order part of the Poisson equation [2.4] to eliminate the 
gravitational potential. 

The equation of motion [2.24] represents the basic model for vortices in this paper. 
In analogous work in other flelds, this model is usually called the quasi-geostrophic (QG) 
approximation (e.g., Ghil & Childress 1987). This model can be generalized in several dif- 
ferent ways (e.g., Petviashvili & Pokhotelov 1992; Sutyrin 1994); for example, additional 
correction terms can be included in the geostrophic balance equation [2.22]. However, in 
this paper, we want to elucidate the basic physical properties of vortices in circumstel- 
lar disks. We thus focus our discussion on the simple model implied by equation [2.24]; 
further generalizations of this model are left for future work. 
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2.5 Summary of Approximations 



The equation of motion [2.24] describes vortices in a circumstellar disk with self- 
gravity included. We now define a quantity T, which is a generalized vorticity per unit 
surface density, 

T = :T7^(V'/ii + 47rGpi) + . (2.25) 

We refer to this quantity T as the generalized vortensity. The first part of equation [2.25] 
represents the vortensity due to a perturbation; the second part represents the vortensity 
of the unperturbed disk flow. With this deflnition, our equation of motion can be written 
in the seemingly simple form 

DT = 0. (2.26) 

This equation thus implies that the total generalized vortensity is advectively conserved. 
Vortensity can, in principle, be freely transferred between the mean flow and the pertur- 
bation, but the total vortensity is conserved. 

Before considering the solutions to the equation of motion [2.26], we summarize the 
approximations we have made to obtain this form. In particular, we have made the local 
approximation and the assumption of quasi-geostrophic balance. 

The local approximation assumes that the vortex size A is small compared to the 
radius ro of the disk (at the vortex center). Since we are working with two dimensional 
thin disks, we also require that the vortex is large compared to the disk scale height H . 
Our treatment thus implies the following ordering of length scales: 

if < A < ro . (2.27) 

The approximation of quasi-geostrophic balance assumes that the most important 
part of the perturbed velocity fleld corresponds to vortical motions about the vortex 
center (the origin of local coordinates in our treatment). To see what this approximation 
implies mathematically, consider the force equation [2.16]. The third, fourth, and flfth 
terms of this equation give us the geostrophic approximation [2.22]. Thus, for consistency, 
we must require the ordering 

— ~^<0t; 2.28 
r A ^ ' A ' ^ 

where (y) is an average perturation velocity and r is the timescale over which the velocity 

changes. This ordering requirement can be satisfled provided that (1) the timescale for 

the vortex to change its velocity proflle is long (r large), (2) the mean perturbation speed 

is subsonic (so that (u)^ < /ii), and (3) the mean vortex rotation rate is slower than the 

local rotation speed of the disk (so that {v) / A < O). 



3. RESULTS FOR POINT VORTICES AND LINEAR VORTICES 

In this section, we consider the simplest type of vortices in circumstellar disks. We 
discuss the limiting case of point vortices in §3.1 and general linear vortices without 
self-gravity in §3.2. 
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3.1 Point Vortices 



The simplest type of vortex is a point vortex, i.e., a vortex with vorticity locahzed 
at a point in space. Vortices of this type are the solutions to the simplest form of our 
equation of motion. For example, if we consider the rotation rate O and the total surface 
density a of the disk to both be constant, then the equation of motion [2.26] simplifies 
to the form 

Dcu = 0. (3.1) 

In other words, in this limiting case, the vorticity itself is advectively conserved. This 
equation of motion has solutions which correspond to point vortices. Mathematically, we 
can write these solutions as 

= r(52(x-Xo), (3.2) 

where <^^(x) is the two-dimensional Dirac Delta function and Xq is the position of the 
vortex. The quantity F, which determines the magnitude of the vorticity, is usually called 
the circulation. 

Once the vorticity solution is specified in the form of equation [3.2], we can determine 
the physical structure of the vortex (Marcus 1993). We find the velocity field in a manner 
completely analogous to finding the magnetic field produced by a line-like current (e.g., 
Jackson 1975). This calculation implies 

v=-^(^, (3.3) 
Zttzu 

where we have introduced local cylindrical coordinates centered on the position of the 
point vortex {zu is the local radial coordinate). The density enhancement is determined 
by the condition [2.22] of geostrophic balance. Notice that positive density enhancements 
correspond to negative values of the circulation F. In other words, the vortex motion 
must be clockwise (in the —(p direction) for a positive density vortex. For this case, the 
enthalpy perturbation takes the form 

Oo |r| 

hi = — Iog[ri7o/ri7] , (3.4) 

vr 

where zuq is an integration constant which determines the physical size of the vortex (see 
below). Finally, we can write the vortex solution in terms of the surface density via 



(7i _ /^o 

(To \ W 



7 

(3.5) 



7 



where the index 7 is defined by 



^oiri , , 

7^^. 3.6 

Notice that there exists a contradiction in the above argument. We have assumed that 
the total surface density is a constant in order to obtain the simplified equation of motion 
[3.1]. On the other hand, our solution implies a surface density perturbation [3.4 - 3.6]. 
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It turns out that this solution is a hmiting case of the more general problem considered 
in the next subsection. This present solution is valid in the limit that the parameter 
A;/j — )• 0, where is the Rossby wavenumber defined below (equation [3.10]). 

3.2 Linear Vortices Without Self- Gravity 

For the next higher order of approximation, we consider the case of vortices with 
no self-gravity and we consider only the leading order perturbations for O and a. We 
also specialize to the case of a purely Keplerian rotation curve, i.e., O ~ r~^/^ and 
rc?xO = —(3/2)0. In this approximation, we introduce a stream function $ defined by 

Next, we separate the total surface density a into the unperturbed part co(r) and the 
perturbation ai . Furthermore, the unperturbed surface density itself can be expanded as 

(Jq(x) = (7o(r) + X = (7o(r)\l - px/r] , (3.8) 

dr 

where we have assumed that the unperturbed surface density has a power-law distribution 
with index p. Recall that r is the radial coordinate of the vortex center. To leading order, 
the generalized vortensity can be written 



T 



1 



V2$ + l(0-Oo)-%$ + Oo^ , (3.9) 
z a/ Zr . 



where we have subtracted out the constant contribution, ~ Oq/cq, of the background 
state. Keep in mind that we have assumed that the perturbed surface density is small, 
i.e., |cri| <^ ctq. In §4.2, we show how to generalize this result. 

Departing slightly from the usual conventions (e.g., see Marcus 1993), we define a 
Rossby wavenumber Ur^ i.e., 

kR = ^. (3.10) 



a 



For circumstellar disks, this Rossby wavenumber is the inverse of the disk scale height 
(kR = 1/H). We can also expand the term O — Oq to leading order in x = r — tq. We 
thus obtain the expression 



T ' 



- kl^ + aox 



(3.11) 



where we have defined a parameter ckq, 

ao = Oo(2p- 3)/4r . (3.12) 

Notice that in general ckq can be either positive or negative, depending on whether the 
surface density co(r) decreases faster or slower than the rotation curve 0(r). For the 
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particular case of p = 3/2, we obtain ckq = 0. Theories of the formation of protosteUar 
disks suggest that disk density profiles will have indices in the range 1 — 2 (see, e.g., 
Cassen & Moosman 1981; Terebey, Shu, & Cassen 1984). Thus, we expect ckq to be small 
and we can take p ^ 3/2 and hence ckq ~ as a starting approximation. 

The solutions for the stream function, the velocity field, and the surface density 
perturbation are qualitatively similar to those of the point vortices discussed in the pre- 
vious subsection. We can find these solutions as follows. We first construct generalized 
vortensity profiles T that satisfy the conservation condition of equation [2.26]. We can 
then invert the differential operator appearing in equation [3.11] to obtain the stream 
function $. Given the stream function, we can determine the velocity field of the vortex 
through the relation 

v = 5xV$. (3.13) 

For example, suppose we take a solution T(x) for the vortensity, which could be a collec- 
tion of point vortices, or perhaps a constant over some area. We thus obtain the solution 
for the stream function in the form 

$(x) = (7o(r) J (/2x'G(x,x')T(x'), (3.14) 

where we have taken ckq = and where G(x,x') is the Green's function corresponding 
to the differential operator in equation [3.11]. This operator is simply the modified 
Helmholtz operator (in two dimensions). The appropriate Green's function is made up 
of modified Bessel functions (e.g., Arfken 1985) and can be written in the form 

G(x, x') = ^ KoikRl^ - x'l) , (3.15) 
zvr 

where A'o is the modified Bessel function of order zero. Note that this is the appropriate 
Green's function for the boundary condition that G — )• as |x — x'| — )• oo. Notice also that 
the Rossby wavenumber determines the length scale on which the stream function $ 
changes; the velocity field v changes on this same length scale. The lengthscale 1/kR is 
comparable to the scale height H in the disk; thus, the ordering requirement of equation 
[2.27] can be satisfied, but only by a small margin. 

Thus, for a given solution T(x) for the vortensity, we can write the solution for the 
stream function in the form 

^/^o = ^ f d'^'Ko{kR\^ - x'l)T(x') . (3.16) 

For the purpose of illustration, let the generalized vortensity be that of a point vortex, 
i.e., 

r = —6^x). (3.17) 

For a brief discussion of more general solutions for the vortensity T, see Appendix A. 
With this choice [3.17] for T, the stream function becomes 

$ = ^Ko(kRw) = (3.18) 
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where vd is the local radial coordinate centered on the vortex. Equation [3.18] thus 
determines the density perturbation. The velocity field of the vortex produced by the 
stream function can be determined from equation [3.13], which implies that 

^ = ^ X VKo(A;Kr^)] = {-KA^Rw)\ (p , (3.19) 

where we have used the relation Kq(z) = —Ki(z) in obtaining the second equality. This 
solution is thus similar to the point vortex solution found in the previous subsection. In 
the limit of small local radius, — )• 0, we find that v oc T/w (compare with equation 
[3.3]). However, in the opposite limit w — )• oo, we find v oc (TkR)z~^/'^e~^ ^ where 
z = kRW. The perturbation is thus exponentially suppressed at large radial distances 
from the vortex center, i.e., for Urw ^ 1. 

Next, we consider vortices which have a constant vortensity over some area (and 
zero vortensity outside). Whenever the shear fiow can be neglected (so that we can 
take O = Oo), we obtain circularly symmetric vortices. Here, we consider vortices with 
vortensity profiles which are step functions in radius, i.e., we consider the form 

T = -|^0(i?-r), (3.20) 



vr 

where is the usual Heaviside step function and where R is the radius of the area with 
nonzero vortensity. Using this form in equation [3.14], we can directly calculate the 
stream function, the density perturbations, and the corresponding velocity fields. The 
results are shown in Figures 3 and 4 for the cases with UrR = 1/5 and 5. The solutions 
for the velocity field and the vorticity can be found analytically (Marcus 1993). For small 
radii, < w < we obtain 

Y = ^Ki(kRR)h(kRw)^, (3.21) 

and 

^ = ^ K,{kRR)Io{kRw) , (3.22) 
ttR 

whereas for larger radii, w > R^ we obtain 

Y= ^ h{kRR)R\{kRm)^, (3.23) 
ttR 

and 

Tkn 

to = ^ h{kRK)K^{kRw) , (3.24) 

vrix 

where 1^ and are modified Bessel functions. In the limit of small kRR, these solutions 
approach the point vortex solutions discussed above. For kRR ^ 1, the vortices have 
"quiet centers" and most of the vorticity is concentrated near the boundary. 

In the presence of shear fiow (neglected here), vortex solutions with constant T 
become noncircular. For the special case of linear shear and kR — )• 0, the vortices have 
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the shape of perfect eUipses (Moore & Saffman 1971). For the more general case, the 
vortices are completely specified by the location of their boundaries, but the boundary 
shapes must be found numerically (for further discussion, e.g., see Marcus 1993). 



4. HIGHER ORDER EFFECTS 

In this section, we consider higher order effects on the vortices. In particular, we 
consider the effects of self-gravity (§4.1), nonlinearity (§4.2), magnetic fields (§4.3), and 
the difference between vortices with positive and negative circulations (§4.4). We also 
briefly address the important issue of vortex generation in §4.5. 

4-1 Leading Order Effects of Self- Gravity 

In order to understand the effects of self-gravity on vortices, we proceed in the same 
fashion as in the previous subsection but we include the gravitational term. Repeating 
the steps leading to equation [3.11], we obtain 



T 



1 







(4.1) 



where we have again taken ckq = 0. We now make a heuristic approximation to simplify 
the gravity term, 

Pi = — = -TT^hi = —TTT^ ' (4-2) 



Z d^Z d^Z 

where Z is an effective total height of the disk. We also deflne an effective Jeans wavenum- 
ber kj through the relation 

2 _ ArrGao 



The expression for vortensity becomes 

T 1 



(To(r) 



dZ 



- + A;^$ 



(4.3) 



(4.4) 



We can now directly see the effects of gravity on the structure of the vortex. In the 
previous example, we saw that the Rossby wavenumber determines the lengthscale on 
which the vortex structure can change. Gravitational effects appear in the equation with 
the opposite sign; as a result, gravity makes the effective Rossby wavenumber smaller 
and hence makes the effective Rossby radius larger. 

An important crossover point occurs when the gravitational contribution exceeds 
the (old) Rossby wavenumber contribution and the total effective Rossby radius becomes 
imaginary. This crossover occurs when = kj, which implies that 



47rGcro 



^Iz . 



(4.5) 
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If we also assume that the total effective height Z of the disk is twice the usual thermal 
scale height in the disk, i.e., Z 2alQ,Q, we obtain the crossover condition in the form 

(5t~2, (4.6) 



where Qt is the Toomre stability parameter for gaseous disks (Toomre 1964). Here, when 
Qt becomes less than 2, the effective Rossby wavenumber becomes imaginary. Notice 
that in identifying our crossover condition with Qt^ we have assumed that the rotation 
curve is purely Keplerian so that the epicyclic frequency k, is equal to the rotational 
frequency O. 

The parameter Qt must be larger than unity in order for the disk to be stable to 
axisymmetric disturbances. In fact, consideration of non-axisymmetric perturbations in 
circumstellar disks suggests that the stability parameter must be (at least) as large as 
Qt ~ 2; otherwise the disks would be highly unstable and would have short lifetimes 
(e.g., Adams, Ruden, & Shu 1989; Shu et al. 1990; Tomley, Cassen, & Steiman-Cameron 
1991; Adams & Benz 1992; Laughlin 1994). As a result, many disk systems tend to live 
close to this crossover point. As disks evolve in time, they eventually lose mass to the 
central star, the value of Qt increases, and the effective Rossby wavenumber approaches 
the true Rossby wavenumber. 

The effect of self-gravity is thus to make the effective Rossby radius longer. It is 
useful to define an effective Rossby wavenumber A^eff, 

A;eff^^4{l-7|-}- (4.7) 
Qt 

The equation of motion which determines the stream function can then be simplified to 
the form 

V'$- A;eff'$ = aoT(x). (4.8) 

Solutions for $ can be found in a manner analogous to that presented in the previous 
section. 

For purposes of illustration, we consider the special case in which the effective Rossby 
wavenumber vanishes. In this particular case, the appropriate Green's function for the 
differential operator in equation [4.8] (i.e., a two dimensional Laplacian operator) is given 
by 

G(x,x') = -:^log|x-x'|. (4.9) 
The solutions for the stream function are thus given by 

$ = -— /(/2x'T(x')log|x-x'|. (4.10) 
2vr J 

Next, we consider the simplest case in which the vortensity is localized and can be 
represented as a delta function as in equation [3.17]. We thus obtain 

$ = -i^log|x|, (4.11) 
zvr 
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where vd is the local radial coordinate and the velocity is in the local azimuthal direction 
(p. Thus, the velocity field of the perturbation is the same as that of the point vortices 
studied in §3. 

1^.2, Nonlinear Density Effects 

In this subsection, we generalize our previous results by considering the case when 
the perturbed surface density ai is not small compared to the unperturbed state. For an 
isothermal equation of state, the enthalpy h is related to the surface density through the 
relation 

/i = aMogcr, (4.13) 
which implies that the perturbed part of the enthalpy can be written in the form 

hi = h - ho = a^log — . (4.14) 

l(7o J 

Using this result, we can write the generalization of equations [3.9] and [3.11] in the form 

+ lo + = T(x)ao exp[20o$/a'] • (4.15) 

Z iLQ 

In order to isolate the nonlinear effects in this equation, we consider the particular 
case in which the surface density ctq and the sound speed a are all constant. This ap- 
proximation is valid when these quantities vary sufficient slowly compared to the stream 
function. We define a reduced stream function according to 

20o hi 

/i = ^$ = 4- 4.16a 
We also define a non-dimensional variable z by 

z = Urx. = — . (4.16b) 
With these definitions, the equation of motion for the vortensity can be written 



2 _ 2Tcro ^ 47rGao 

Oo " O^Z ^ ^' (4.17) 



=Ae'' - Bfi-1. 



The quantity A defined above is the ratio of the vortensity in the vortex to that of 
the unperturbed circumstellar fiow. The second constant B is the ratio of the Jeans 
wavenumber (defined in equation [4.3]) to the Rossby wavenumber (defined in equation 
[4.4]); thus, the constant B determines the relative importance of self-gravity. 
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This nonlinear form of the equation of motion can lead to qualitatively new behavior 
under some circumstances. In the general case (see Appendix A), the vortensity T is some 
function of the stream function ($ or /i), so that in equation [4.17], the quantity A = 
Thus, the right hand side of equation [4.17] can be considered as some function F(fj,). In 
other words, the equation becomes 

dV 

V^/i = A(/i)e'' - 5/i - 1 = = - — . (4.18) 

The function -F(/i) thus plays the role of a force and its integral V(fi) plays the role of 
a potential. Nonlinear equations of this type are well studied in the context of quan- 
tum field theory (e.g., Coleman 1985; Rajaraman 1987). The properties of the function 
V(fi) determine the allowed types of solutions for /i(z). In particular, when V has both 
a double zero and a simple zero (when considered as a function of the reduced field 
/i), then solitary wave solutions are allowed (see, e.g., Drazin & Johnson 1989; Infeld 
& Rowlands 1990); these solutions correspond to spatially localized wavelike structures. 
Unfortunately, however, a complete analysis of all of the possible solutions to the nonlin- 
ear equation of motion is beyond the scope of this present work and must be left for the 
future. 

4-3 Effects of Magnetic Fields 

In this section, we consider the effects of magnetic fields on these vortices. For the 
sake of definiteness, we will consider a magnetic field which threads the disk in the vertical 
direction only, i.e., 

B = B{x,y)z. (4.19) 

We note that the no-monopole condition, V • B = 0, is automatically satisfied with this 
type of magnetic field. The field adds a force term to the equations of motion; this force 
fjw has the form 

fM = ^(V X B) X B = -^VB^ . (4.20) 

Thus, one important effect of the magnetic field is to add an additional contribution to 
the pressure force. For the particular field geometry used here, the magnetic field simply 
changes the effective equation of state of the fluid. The change is thus quantitative, but 
not qualitative. 

For completeness, we note that the continuity equation for the magnetic fleld can be 
written in the form 

B 



D 



0. (4.21) 



In other words, the magnetic fleld strength per unit surface density is also advectively 
conserved. 

The presence of a magnetic fleld also allows for a new channel of energy dissipation. 
In general, if the fluid has a flnite electrical resistivity the magnetic fleld will tend to 
diffuse according to the equation of motion 

— + V X (B X u) = r/V^B. (4.22) 
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In the limit of — )• 0, we recover the usual flux freezing assumption. For a given resistivity 
?7, the magnetic diffusion timescale tm is given by 

tm-AV??, (4.23) 

where A is the physical lengthscale of the system (i.e., the vortex size). Using typical 
values for rj (see Shu 1992), we flnd that the diffusion timescale is much longer than the 
expected disk lifetime (~ 10^ — 10^ yr, e.g., Strom, Edwards, & Skrutskie 1993). 

4-4 The Prograde/Retrograde Asymmetry 

Within Keplerian circumstellar disks, vortices with positive circulation behave dif- 
ferently than vortices with negative circulation. As we have discussed previously, in the 
present context, a vortex with a positive density perturbation must have a negative cir- 
culation in order to achieve geostrophic balance. In other words, a vortex with a positive 
density perturbation must rotate clockwise in the rotating frame of reference in which 
we have been working. Now, suppose a planet (or other secondary body) forms within a 
vortex. This planet would naively appear to be rotating in a retrograde sense. However, 
this apparent retrograde rotation (say, at a rate Op) must be corrected for the fact that 
the frame of reference is rotating as well in the opposite direction (at the local Keplerian 
rotation rate O). As long as the "local" rotation rate Op is smaller than the Keplerian 
rotation rate O, the secondary body will rotate in a prograde sense in the inertial frame 
of reference. Notice also that if Op > O (which would lead to a real retrograde rotation 
of the planet), then the approximations leading to our equations of motion break down. 

Another important difference between vortices with positive and negative circulation 
is their lifetime and stability. In a Keplerian disk, the mean flow velocity decreases with 
radial distance to the star. In our local frame of reference centered on the vortex, the 
mean flow velocity is negative for x > and positive for x < 0; this effect is best illustrated 
pictorially (see Figure 5). As a result, vortices with negative (clockwise) circulation rotate 
"with the flow", whereas vortices with positive (anti-clockwise) circulation must directly 
oppose the flow; these latter vortices are deflned as "adverse". Detailed numerical studies 
(cf. the review of Marcus 1993) show that adverse vortices are quickly ripped apart into 
long fllaments and fragments, i.e., their lifetimes are very short. On the other hand, 
vortices of the other type tend to merge and grow into coherent structures. A complete 
discussion of the lifetime and stability of vortices is beyond the scope of this present 
paper. We stress, however, that the vortices with positive density perturbations (negative 
circulation) are the ones most likely to live a long time in a Keplerian flow. 

4-5 Vortex Generation 

Perhaps the most important unresolved issue is the manner in which these vortices 
are generated in the flrst place. Thus far in this paper, we have found many different 
solutions to the equations of fluid dynamics; these solutions represent possible behavior of 
the circumstellar disk. However, before the astrophysical importance of vortices can be 
established, one must understand the generation mechanism. In particular, the timescales 
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for vortices to grow and decay are important. In this section, we (crudely) estimate a 
timescale for the vortensity (or vorticity) to change in a circumsteUar disk. 

One can argue that aU differentiaUy rotating fluid systems naturaUy produce vortical 
motions and hence vortices must be important at some level. For example, the Earth's 
atmosphere, the Earth's oceans, and the atmospheres of the giant planets all produce 
many different kinds of vortices such as those studied here. In addition, many types of 
vortices can be easily generated in laboratory experiments (Nezlin & Snezhkin 1993). 
Thus, one might naively expect vortices to arise naturally in circumsteUar disks. 

One important constraint on vortex generation is provided by Kelvin's circulation 
theorem (e.g.. Landau & Lifshitz 1987; Ghil & Childress 1987; Shu 1992). This theorem 
states that for inviscid barotropic flow^ the number of vortex lines that thread a given 
area (that moves with the fluid) remains unchanged with time. Thus, for the case of 
circumsteUar disks which are nearly barotropic, vorticity is not automatically "generated" 
in these systems. In the context of our local equation of motion (see equation [2.20]), the 
total vorticity (per unit surface density) is advectively conserved. This total includes both 
the vorticity of the perturbation and that of the mean flow. Thus, vorticity can be freely 
exchanged between the mean flow and the perturbation, but cannot be directly generated 
(or destroyed). We also note that, except under rather special circumstances, the mean 
flow is not unstable and does not spontaneously transfer its vorticity into perturbations 



If the flow is not exactly barotropic, then vorticity can be generated directly. Con- 
sider, for example, the force equation [2.16] with a standard pressure term of the form 
p~^Vp. When we take the curl of the force equation to obtain the equation of motion 
for the vorticity, we obtain a forcing term of the form 



where we have used the ideal gas law p = pRT in obtaining the second equality. Clearly, 
for a barotropic equation of state p = p{p), this new term vanishes. However, this forcing 
term will be nonzero whenever the surfaces of constant temperature do not line up exactly 
with the surfaces of constant density. Such a situation can occur in a circumsteUar disk. 
In the absence of disk accretion energy, one important heating source for the disk will be 
reprocessing of stellar photons; this energy source naturally leads to nearly axisymmetric 
surfaces of constant temperature. On the other hand, gravitational instabilities naturally 
produce spiral patterns and hence non- axisymmetric surfaces of constant density (e.g., 
Adams, Ruden, & Shu 1989; Laughlin & Bodenheimer 1994). 

We now consider the idealized case in which the surfaces of constant temperature 
are exactly axisymmetric. We will also assume that the disk has strongly growing spiral 
density perturbations so that the surfaces of constant density are tilted at an angle i with 
respect to axisymmetric surfaces. Spiral density wave theory in the WKBJ limit (see Shu 
1992; Binney & Tremaine 1987) implies that the tilt angle is given by 



(vortices). 



Fy = — Vp X Vp = — Vp X VT , 
P P 



(4.24) 



m 

tan I ~ - — , 

kr 



(4.25) 
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where m is the azimuthal wavenumber of the perturbation and k is the radial wavenumber. 
For global modes in circumstellar disks, we expect those with m = 1 to be the most 
important (Adams, Ruden, & Shu 1989; Laughlin & Bodenheimer 1994). We also expect 
the radial wavelength to be comparable to the radius and hence i ~ l/27r. Given this set 
of approximations, the generalized equation of motion for the vortensity can be written 

DT ^ -^-\Vp\ |VT| - Fdis , (4.26) 
Zttct p 

where we have taken sinz ~ z ~ l/27r. For completeness, we have also included a 
dissipation term Fois- 

Next, we want to estimate the timescale Ty for vortensity in the disk to be generated 
by the effect described above. To make this estimate, we write DT ~ co/aTv- We write 
the density gradient term as \\^p\/p ~ where /3 is the amplitude and L is the 

sizescale of the spiral perturbation. Finally, we write the temperature gradient term as 
|VT| ~ qT/r^ where q is the power-law index of the temperature profile. Solving for the 
timescale, we obtain the estimate 

27r LorL 

• (*-^') 

where a is the sound speed. Very roughly, we expect /3g ~ 1, coA ~ a, and L ~ A, so 
that the timescale for the vortensity to change is Ty ~ 27rr/a ~ 500 yr. This timescale 
corresponds to several orbit times and is thus interesting for the evolution of circumstellar 
disks. We stress, however, that this argument is extremely crude and must be improved 
with an honest calculation. 

Before leaving this section, we note that several recent papers have begun to address 
the issue of vortex formation in disks and other rotating systems. For example, the 
generation of large-scale flow patterns in turbulent astrophysical systems has been studied 
by Kitchatinov, Riidiger, & Khomenko (1994); they flnd that when the rotation rate 
becomes sufficiently large, the system becomes unstable to the production of large-scale 
motions which can be considered as vortices. Related work has studied the emergence of 
coherent vortices out of turbulent flows (McWilliams 1984, 1990). 



5. ASTROPHYSICAL APPLICATIONS 

In this section, we discuss several different astrophysical applications of vortices. The 
most important applications are in circumstellar disks, where vortices may enhance the 
formation of giant planets (§5.1) and may provide enhanced energy dissipation for disk 
accretion (§5.2). For the sake of completeness, we also briefly consider applications of 
these vortices to the formation of molecular clouds out of the galactic disk (§5.3) and to 
the formation of substructure within molecular clouds (§5.4). 

5.1 Dust Separation and the Formation of Giant Planets 
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One of the main difficulties in forming giant planets (Jupiter, Saturn, Uranus, Nep- 
tune) through the mechanism of gravitational instability is that the planets are highly en- 
riched in heavy elements (see, e.g., DeCampli & Cameron 1979; Stevenson 1982; Cameron 
1988; Podolak, Hubbard, & Pollack 1993). In the solar nebula from which the planets 
formed, most of the heavy elements are expected to be in the form of dust grains. For 
most planet formation scenarios considered previously, the timescale for the dust grains 
to migrate to the center of the forming protoplanet is too long to explain the observed 
enrichment. In addition, a solid rocky core is found at the center of the giant planets; 
it is difficult for such a core to form at all in most previously considered models. As we 
show below, however, vortices can, in principle, provide a mechanism for separating dust 
and gas on a shorter timescale and may provide a means of forming giant planets within 
the gravitational instability scenario. 

In a rotating vortex system, the dust grains are not completely coupled to the gas. 
The grains do not directly feel the pressure forces, but they are coupled indirectly through 
drag forces. Following convention, we define a response time Td which determines the 
timescale on which dust grains respond to the force exerted by gas drag. For the regime 
of parameter space appropriate for circumstellar disks, the response time is given by 

Td = - , (5.1) 

Zpa 

where (h) is the average radius of the dust grains and pd is the grain mass density (see, 
e.g., Weidenschilling & Cuzzi 1993). For typical conditions {pd ~ 2 g/cm^, (h) ~ 10""^ 
cm, p ~ 10""^^ g/cm^, a ~ 5 X 10"^ cm/s), this timescale is relatively short, Td ~ 2000 
s. Under these conditions, the dust grains quickly reach their terminal velocity, i.e., the 
force due to gas drag is balanced by the inward radial force (here due to Coriolis effects). 
This terminal velocity Uterm is given by 

dhi 

t^term = ^d^— • (5.2) 

aw 

In order to estimate the timescale for dust grains to settle to the center of the 
vortex, we use the point vortex solution found analytically in §3. If we use this solution 
in equation [5.2] and solve the resulting differential equation, we obtain the following 
simple expression for the timescale Tgettle, 

_ TT 

2r,0o|r| ' ^^-^^ 

where A is the assumed starting radius of the dust grain and is taken to be the vortex 
size. To obtain a numerical estimate for the timescale, we let |r| ~ aA (corresponding 
to a weakly nonlinear vortex) and we let the vortex size be the inverse of the Rossby 
wavenumber, i.e., A ~ a/Oo- We thus obtain the dust settling timescale 

Tsettie = 2 X lOVr (Oo lOOyr)-^ . (5.4) 
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This timescale is comparable to that expected for dust setthng in a hydrostaticaUy sup- 
ported protoplanet. Furthermore, this timescale is too long to explain the observed 
element segregation in the giant planets in our solar system. In order to overcome this 
difficulty, either the grain size must be larger or the vortex strength (as measured by the 
magnitude of the circulation |r|) must be larger than we assumed here. For example, 
suppose we require the dust settling timescale to be ~ 10^ yr, roughly two orders of 
magnitude shorter than given above. This time requirement can be met provided that 
mean dust radius (6) and the circulation F satisfy the following constraint: 

(b) iri , , 

f^^>100. 5.5 

i/im a,/kR 

Since we naively expect both of the ratios in the above equation to be of order unity, 
we find that sufficient dust segregation to form giant planets requires somewhat extreme 
conditions. We note, however, that in principle both the circulation and the grain size 
can be large and hence giant planet formation is at least possible is this scenario. 

The above simple calculation can be generalized to obtain a more accurate descrip- 
tion of dust settling in vortices. Recent papers (Tanga et al. 1995; Barge & Sommeria 
1995) have performed numerical integrations of the equations of motion for dust particles 
in vortices and confirm that dust particles concentrate inside vortices on a relatively short 
timescale. 

Finally, we note that vortices enhance the possibility of forming giant planets via 
gravitational instability by making the local surface density larger. The perturbation 
requirements to form a secondary body in a circumstellar disk can be written in fairly 
general form; this argument is given in Appendix B. The basic result is that a moderately 
large perturbation, with density contrast /3 = Aa / a ~3-5, is required in order to form a 
secondary body. 

5.2 Energy Dissipation and Disk Accretion 

In this subsection, we briefly discuss how vortices can affect energy dissipation and 
hence accretion in circumstellar disks. As we discuss below, the overall effect of vortices is 
to provide additional avenues for energy dissipation in disks. We can identify three con- 
ceptually different ways for vortices to participate in energy dissipation: (1) annihilation 
of vortices, (2) energy dissipation within a single vortex, and (3) a cascade of vortices. 

The annihilation of vortices directly leads to the dissipation of energy. The simplest 
model to consider is that of point vortices with circulation of opposite sign (see §3). The 
vortices can annihilate, or "short out", in much the same way that magnetic fleld lines 
of opposite direction dissipate their energy. 

In order to consider the dissipation of energy within a single vortex, we think of the 
vortex as a small scale analog of the circumstellar disk itself. As shown from our solutions 
to the equations of motion in the previous sections, the vortices are strongly differentially 
rotating. Now suppose the fluid has a viscosity i^, which we parameterize in the usual 
way according to 

2 

ly = —aaH ~ aaA , (5.6) 
3 
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where A is the characteristic vortex size and a is the usual dimensionless "a parameter" 
(see, e.g., Pringle 1981). The viscous diffusion timescale td is then given by 



A' A/« 1 6yr 
V a ila a 

where we have taken typical values to obtain the numerical estimate. Many studies of 
viscous disks suggest that a ~ 10~^ — 10""^ (e-g-? Lin & Papaloizou 1980, 1985; Lin 1981), 
which implies a vortex dissipation timescale of rj^ ~ 6 X 10^ — 6 X 10"^ yr. 

Because of the smaller size scale of the vortex, the energy dissipation rate is enhanced 
over the accretion timescale of the disk as a whole by a factor / given by 

/ ~ (A/r)^ ~ 100. (5.8) 

Thus, for a given fluid viscosity, vortices are more eflicient at dissipating energy than the 
disk as a whole. As a result, vortices can play an important role in energy dissipation and 
hence disk accretion if a mechanism exists to efliciently transfer energy from the mean 
flow into the vortex (see §4.5). 

Finally, we consider the possibility that these vortices could represent the flrst stage 
in a cascade of eddies. We have shown that vortices can arise from a circumstellar 
disk. We have also shown that the vortices themselves are strongly diff'erentially rotating 
systems. As a result, smaller vortices can, in principle, arise within the diff'erentially 
rotating flow of the the larger vortices. These smaller vortices can support still smaller 
vortices and so forth. In this manner, a cascade of eddies can be set up and can transfer 
energy from larger scales down to smaller scales. This energy transfer is exactly what is 
needed to provide an anomalous viscosity in the disk. We stress, however, that we still 
do not understand the mechanism(s) by which vortices are generated; as a result, this 
cascade of eddies remains little more than a speculative possibility. 



5.3 Molecular Cloud Formation from the Galactic Disk 

The disk of the galaxy is in many ways a larger version of the circumstellar disks 
that we have been considering thus far. In particular, vortex motions can arise in the 
galactic disk and can form substructure within the galaxy (see Petviashvili & Pokhotelov 
1992 for a review). The treatment is analogous to that presented above. Although the 
rotation curve of the galaxy is non-Keplerian, to leading order, only the local rotation 
speed is important. In any case, the relevant size scale for any structures produced by 
vortices in the galaxy will be given by the Rossby radius A = cr*/0, where O is the local 
rotation speed of the galactic flow and where cr* is the velocity dispersion. For typical 
values for our galaxy, O ~ 10""^^ rad/s and cr* ~ 30 km/s (e.g., Binney & Tremaine 
1987), we obtain 

^^^^-^ = ' (soSv^) (lo-w/s) ' • ^^-^^ 

Notice that the flnal size of the structure can be much smaller than this value (after some 
collapse has taken place). The corresponding mass scale is also fairly large, Mc ~ SA^ 
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~ 10 — 10^ However, if the cloud formation process is not completely efficient and 

only about 10% of the mass ends up in the cloud, then this mass scale is about right for 
giant molecular clouds. 



5.^ Formation of Molecular Cloud Substructure 



It is possible for vortical motions to help produce substructure on the size scale of 
molecular clouds. The two-dimensional vortices studied in this paper are not well suited 
to this problem, as molecular clouds are inherently three-dimensional. We note, however, 
that any type of rotating fluid system supports vortical motions and hence vortices can 
occur in these clouds. Very roughly, the relevant size scale for any possible substructure 
is determined by the Rossby radius A = a/0. Using typical values for the temperature 
and rotation rates of molecular clouds (e.g.. Blitz 1993), we flnd 



This size scale is reasonable for molecular cloud substructure. A similar type of self- 
gravitating vortex has been considered in the context of forming substructure in the 
interstellar medium (Chantry, Grappin, & Leorat 1993; see also Sasao 1973). In addition, 
the importance of vorticity in the collapse of molecular clouds and the formation of 
extended fllamentary structures has recently been studied by Monaghan (1994). 



In this paper, we have begun a study of vortices in circumstellar disks. These vortices 
are basically storm systems in a generalized geostrophic balance, i.e., a balance between 
pressure forces, gravitational forces, and the Coriolis force. 



We have obtained a number of basic results concerning the physical properties of 
vortices in circumstellar disks and other astrophysical systems. 

[1] Many different types of vortex solutions are possible in circumstellar disks (see §3-5 
and Appendix A). In other words, circumstellar disks can exhibit many different 
types of vortical motions. 

[2] Point vortices are the simplest type of vortex. Their properties can be found ana- 
lytically (see §3) and hence these vortices provide a prototype for understanding the 
basic physics of these systems (see also the review of Marcus 1993). 

[3] In the absence of self-gravity, the relevant sizescale A for this type of vortex is 
roughly determined by the inverse of the Rossby wavenumber, i.e., A = a/0. For 
circumstellar disks, this sizescale is comparable to the thermal scale height H in the 
disk. 




(5.10) 



6. SUMMARY AND DISCUSSION 



6.1 Overview of Results 
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The leading order effect of self-gravity is to make the effective Rossby wavenumber 
smaller (§4.1). This effect, in turn, makes the size of the vortex larger than in the 
case without gravity. 

Nonlinear effects can lead to qualitatively new behavior (§4.2). In particular, under 
special circumstances, solitary wave solutions can arise. 

Magnetic fields can, in principle, have two different effects on the evolution of vortices 
(§4.3). The fields exert a force on the fluid and hence add an additional force term 
to the equation of motion. For the simplest case of a magnetic fleld in the vertical 
(z) direction, this force corresponds to a magnetic pressure only. The second effect 
is to allow for ohmic dissipation; however, for expected parameters in circumstellar 
disks, this effect is small. 

Circumstellar disks can, in principle, generate vorticity (or vortensity) through baro- 
clinic effects (§4.5). These effects arise whenever the surfaces of constant density do 
not line up with the surfaces of constant temperature. Disks can naturally obtain 
nearly axisymmetric temperature proflles (from reprocessing of stellar photons) and 
non-axisymmetric density proflles (from self-gravitating spiral modes). The esti- 
mated timescale for the vortensity to change is ~500 yr. 

Vortices can enhance the formation of giant planets through the mechanism of grav- 
itational instability (§5.1). In particular, dust grains settle to the center of these 
vortices. In order for the timescale for dust settling to be interesting for planet for- 
mation, either the dust grains must be quite large ((6) ^ 1 /xm) and/or the vortex 
circulation must be very strong (|r| ^ aH). Vortices also enhance the formation of 
giant planets by gravitational instability by increasing the local surface density of 
the disk. 

We have derived a general argument which illustrates the conditions under which 
secondary bodies can form in circumstellar disks (Appendix B). For moderately 
massive disks associated with young stellar objects, a moderately large perturbation 
amplitude /3 = Aa / a ~ 3-5 is required to form a secondary. 

Planets formed in vortices should generally rotate in a prograde sense with respect 
to their orbits (§4.4). In order to produce a planet with retrograde rotation, the 
circulation of the vortex must be sufflciently strong that the approximations used in 
this paper break down. 

Vortices can enhance the dissipation of energy in circumstellar disks and can thereby 
help disk accretion (§5.2). For a given fluid viscosity, vortices dissipate energy faster 
than the disk as a whole. The physical reason for this enhancement is that the dif- 
ferential rotation in the vortex is strong and the size scale is short compared to the 
entire disk. Thus, vortices can enhance disk accretion, provided that some mecha- 
nism exists to transfer energy from the mean circumstellar flow into the vortices. 

These vortices can, in principle, also form structure in the galactic disk and within 
molecular clouds themselves (§5.3 and §5.4). The size of these structures is deter- 
mined by the local Rossby radius. In the case of the galactic disk, the mass scale of 
the expected structures is roughly an order of magnitude larger than that of giant 
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molecular clouds. In the case of rotating molecular clouds, the length scale is ~1 pc, 
an interesting and relevant sizescale for substructure in molecular clouds. 

6.2 Disk Accretion vs Planet Formation 

In this paper, we have suggested that vortices may be important for both planet 
formation and for disk accretion. An analogous dual role has been claimed previously 
for gravitational instabilities in disks (e.g., Adams, Ruden, & Shu 1989), although the 
secondary bodies are more often considered to be "binary companions" for this latter 
mechanism. In either case, it might seem paradoxical that one physical process can 
produce two very different results - namely secondary bodies and disk accretion. However, 
this apparent contradiction can be resolved for both mechanisms as follows: Energy 
dissipation and hence disk accretion arises for "failed" structures, whereas secondary 
bodies can form only within "successful" structures. 

We first consider energy dissipation and disk accretion. Self-gravitating spiral modes 
lead to disk accretion when they saturate at a fairly low amplitudes (Laughlin & Boden- 
heimer 1994), i.e., when the instabilities "fail" to achieve strong nonlinear growth. Sim- 
ilarly, vortices can lead to energy dissipation and hence help disk accretion when they 
dissipate their energy faster than they can grow, i.e., when they also "fail" to achieve 
strong nonlinear growth. 

On the other hand, secondary bodies can be produced by either mechanism for 
the case of "successful" structures. Self-gravitating spiral modes can collapse to form 
secondary bodies when they successfully grow well into the nonlinear regime; as shown 
in Appendix B, a density contrast /3 ~ 3 - 5 is required (see also Adams & Benz 1992; 
Laughlin, Bodenheimer, & Yorke 1995). Vortices can lead to giant planet formation if 
they are sufficiently long-lived (so that dust grains collect at the vortex center) and/or 
highly nonlinear (so that they must collapse gravitationally), i.e., if the vortices are 
"successful" . 

6.3 Future Work 

This paper represents a preliminary step toward understanding the physics of vortices 
in circumstellar disks and related astrophysical systems. As we outline below, many 
directions of future work remain. 

The most important issue to be studied is the manner in which vortices are generated 
in circumstellar disks (see §4.5). This process must be understood in order to determine 
what types of vortices are produced, how often they are produced, and where in the disk 
they are produced. 

Another important issue is the question of vortex stability and lifetimes. In order for 
vortices to be interesting for giant planet formation, they must live for many dynamical 
timescales and hence be relatively stable. On the other hand, vortices which dissipate 
their energy rapidly will affect disk accretion. It is thus important to determine the 
conditions under which vortices are long-lived and short-lived. 

Many higher order issues also remain unresolved. We have considered only a small 
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fraction of the possible vortex solutions. Vortices can exhibit a wide range of complicated 
behavior, including vortex merging, vortex scattering, growth and collapse, decay and 
dispersive spreading, etc. These issues should also be studied. 

Finally, we note that the overall goal of this work is to understand the dynamics 
and evolution of circumstellar disks as a whole. Vortices will play only a partial role in 
determining the long-term fate of disks. One challenge for the future is thus to integrate 
vortices into the overall picture of disk evolution. 
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APPENDIX A: GENERAL SOLUTIONS 
TO THE EQUATION OF MOTION 

We now consider more general solutions to the equation of motion 

DT = 0, (Al) 

which implies that the generalized vortensity is advectively conserved. We specialize to 
the case of steady flow (so that time derivatives are zero) and work in the rest frame of 
the vortex center; we also neglect the background velocity term r(0 — Oq) in the deflnition 
[2.21] of the operator D. With these approximations, the equation of motion becomes 
simply 

v-VT = 0, (A2) 

i.e., the velocity of the vortex flow is perpendicular to the gradient of the vortensity T. 
Next, we note that the velocity is given hy v = z X V$, so that the velocity is also 
perpendicular to the gradient of the stream function $. These two results imply that the 
gradients of the vortensity T and the stream function $ must be parallel, i.e., 

VT = 5V$, (A3) 

where g is an unspecifled function. This equation can be satisfled provided that the 
vortensity is a function of the stream function, 

T = g{<^). (A4) 

To summarize, whenever the generalized vortensity is a function of the stream function 
(equation [A4]), the equation of motion [Al] is satisfled for conditions of geostrophic flow. 
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Thus, instead of choosing solutions for the vortensity T, we can choose a function 
and solve for the stream function from the equation of motion 

V^$-A;|$ + ax = (7oT = (7o^($), (A5) 

which is the appropriate generalization of equation [3.11]. 

We can generalize the above result further by including the background velocity term 
Vb = — ^o)- With this term, the equation of motion becomes 

(1^B^ + v-V)T = 0, (A6) 

where we are again considering the time independent problem. We next define a gener- 
alized stream function $ according to 



$ = $ + y VBdx. (A7) 
The same arguments given above then lead to the generalized solution in the form 

T = g{^) = g{^ + J VBdx), (A8) 
where Q is again an unspecified function. 

APPENDIX B: GENERAL CLUMP ARGUMENT 

We now derive a general argument for the size and mass of clumps that form in a 
self-gravitating circumstellar disk. Suppose that we form a clump in a vortex (such as 
those considered in this paper) or in the collapsing spiral arm of a disk (as in Adams & 
Benz 1992). Let the clump mass be Mc, the clump size be A, and let r be the radial 
location of the clump center. This clump must satisfy three general conditions, which we 
give below. 

(1) The clump size must be smaller than the Roche radius of the star/clump system 
when considered as a binary pair. This condition can be written in the form 

A < ' ^ 



(^)r, (Bl) 



where is the mass of the star. This expression for the Roche radius is valid in the 
limit of small clump mass Mc ^ M*; we have also neglected the gravitational potential 
of the disk itself. 

(2) The clump must be gravitationally bound, i.e., the magnitude of the gravitational 
binding energy must be greater than the thermal energy. Using standard virial theorem 
results (e.g., Shu 1992), we can write this constraint in the form 

A</^?^, (B2, 

a/ 
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where we have used an isothermal equation of state and where the quantity / is a dimen- 
sionless number. For a clump with a modest central concentration, we expect / ~ 1/10. 

(3) Finally, we require that the perturbation which produced the clump was produced 
by the background circumstellar disk. This requirement can be written in terms of the 
amplitude of the perturbation and takes the form 

Mc = /3ao(r)7rA2 , (B3) 

where the parameter jS = a / ao measures the amplitude of the perturbation. 

We can combine the above three constraints to obtain the allowed range of clump 
sizes or masses forming in a circumstellar disk. We find that the allowed mass range is 
given by 

- — > Mr > . (B4) 

As expected, for some range of parameters, the above range of allowed clump masses 
ceases to exist. Also as expected, the requirement for real clump sizes can be written in 
terms of the Toomre parameter Qt as follows: 

Qt < KJI^f' ^ , (B5) 

where n is the epicyclic frequency and is rotation rate due to the stellar potential only, 
i.e., = GM*/r^. For the expected rotation curves in circumstellar disks, k ~ 20*. 
Using representative values for / and we can write our result in terms of the 

required perturbation amplitude for clump formation: 

P>\Qt. (B6) 

Thus, we see that fairly large perturbation amplitudes are required for clump formation. 

This finding is consistent with the results obtained from numerical simulations of 
self-gravitating disks (e.g., Laughlin & Bodenheimer 1994). These simulations show that 
when disks are near stability, gravitational instabilities tend to saturate at somewhat 
low amplitudes (perturbations of order unity) and disk accretion occurs. On the other 
hand, for disks which are more unstable, the gravitational instabilities grow to larger 
amplitudes and the fragmentation of the disk into clumps can occur (see also Bonnell & 
Bate 1994). 
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FIGURE CAPTIONS 



Figure 1. Schematic diagram of local coordinate system. 

Figure 2. Schematic diagram illustrating geostrophic balance. For the case shown, 
both the coriolis force (C) and the gravitational force (g) point radially inward, whereas 
the pressure force (P) is directed radially outward. 

Figure 3. Surface density perturbations for linear vortices (surface density is given 
in arbitrary units). For these solutions, the vortensity profiles are step functions in local 
radius. The radial size of the area of nonzero vortensity is R. Profiles are shown for UrR 
= 1/5 and 5. 

Figure 4. Velocity fields for linear vortices with step function vortensity profiles. 
Notice that the flow is clockwise in the rotating frame of reference moving with the 
vortex center. Cases are shown for (a) UrR =1/5 and for (b) UrR = 5. 

Figure 5. Diagram showing difference between adverse and prograde vortices. When 
the vortical motion of the flow opposes the background shear flow (a), the adverse vortex 
gets ripped apart and has a short lifetime. In the opposite case (b) in which the vortical 
motion goes with the background shear flow, the prograde vortex can be long-lived. 
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